Critical Velocity of Vortex Nucleation in Rotating Superfluid 3 He-A 
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We have measured the critical velocity v c at which 3 He- A in a rotating cylinder becomes unstable 
against the formation of quantized vortex lines with continuous (singularity-free) core structure. 
We find that v c is distributed between a maximum and minimum limit, which we ascribe to a 
dependence on the texture of the orbital angular momentum l(r) in the cylinder. Slow cool down 
through T c in rotation yields l(r) textures for which the measured v c 's are in good agreement with 
the calculated instability of the expected 1 texture. 

PACS numbers: 67.57.Fg, 47.32.-y, 05.70Fh 



A first order transition from one phase to another is 
associated with hysteresis because of the difficulty of nu- 
cleating the new phase. Two effects generally reduce 
the hysteresis. Firstly, thermal or quantum fluctuations 
cause the new phase to appear before the energy barrier 
separating the two energy minima vanishes. Secondly, 
surfaces, impurities, or other external agents reduce the 
energy barrier from its intrinsic value. Both phenomena 
are of crucial importance for the long standing problem of 
critical velocities and vortex nucleation in superfluids jp , 
but occur also in more usual phenomena like formation 
of water droplets or gas bubbles @. The purpose of the 
present work is to study an exceptional case of vortex nu- 
cleation where neither fluctuations nor external surfaces 
should play a role: superfluid 3 He-A. 

In usual superfluids and superconductors the phase slip 
takes place by creation and motion of zeros in the order 
parameter ||. The A phase of 3 He is exceptional be- 
cause the phase slip arises from the motion of the local 
angular momentum axis l(r). The characteristic length 
scale of the l(r) texture is macroscopic ~ 10 /im. There- 
fore all thermal and quantum fluctuations are negligible. 
Moreover, a rigid boundary condition fixes 1 perpendic- 
ular to the wall of the experimental container. Thus the 
processes responsible for the critical velocity take place 
further than 10 /im from the wall, beyond the reach of 
surface roughness. Instead, the critical velocity v c for the 
phase slip depends on the initial 1 texture. We measure 
the critical velocity in a rotating cylinder, and find that it 
may vary within a factor of 6. However, by cooling slowly 
through the superfluid transition temperature T c in rota- 
tion, the equilibrium texture is created and the measured 
v c is in agreement with theoretical calculations. 

Anisotropic superflow. — In ordinary superconductors 
and superfluids, the order parameter has a phase fac- 
tor exp[i(p(r)], and the superfluid velocity is defined as 
the gradient of the phase, v s oc V</>. In 3 He-A there is 
an additional phase factor exp[i0i(p)], which depends on 
the azimuthal angle 4n of the quasiparticle momentum p 
with respect to the angular momentum axis 1. Instead 



of resolving the two phases separately, one may only de- 
fine the total phase factor, which can be expressed as 
(m + in) • p. Here 1, m, and n form an orthonormal 
triad, which generally depends on the location r. The 
superfluid velocity is defined as v s — ^ "^^"-fc' 
where the prefactor equals Planck's constant divided 
by twice the mass of a 3 He atom. This leads to several 
unusual features. For example, let us take an initially 
uniform 1 = z, and then tilt the triads so that l(z) forms 
a helix with an opening angle (3 and a wave vector q. 
This leads to a change in v s . z by ^(1 — cos/3)g without 
a change in the externally applied phase difference [Q. 
Thus 3 He-A can respond to flow by forming an 1 texture. 

The energetics of the current carrying states is con- 
tained in the energy functional [^) 

- Cw ■ v x i + c (i • w)(i • v x i) 

+ iiv s (v • i) 2 + iKt(i • v x i) 2 + \K h \\ x (v x i)| 2 
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(i) 



Here the first two terms describe the anisotropic kinetic 
energy arising from the countcrflow w = v n — v s . The 
terms with coefficients C and Cq contribute to the cou- 
pling between w and inhomogeneous l(r). The five terms 
with Ki coefficients are gradient energies for l(r) and the 
spin anisotropy vector d(r). The last two terms arise 
from the magnetic dipole-dipole interaction, and from the 
external magnetic field H. In the absence of the dipole 
coupling, all uniform current carrying states would be 
unstable [Q. Thus the dipole interaction determines the 
scale of the critical velocity v& = \Jg<i/ P\\ ~ 1 mm/s and 
the length scale £d = 2m \/P||/.9d ~ 10 /im. Except for a 
tiny region near T c (Fig. 2), is much smaller than is 
needed to nucleate usual "singular vortices" , such as one 
encounters in other superfluids 

Let us drive the current by applying a normal fluid 
velocity v n . (This is equivalent to applying a phase dif- 
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ference A(f> = ^v n L between two points a distance L 
apart.) The flow properties depend on the magnitude 
and orientation of the magnetic field - 1 1 . In general, 
a uniform state of 1 is stable for velocities smaller than 
a first critical velocity v c \. At larger velocities there of- 
ten is a stable helical texture. The opening angle /3 of 
the helix grows continuously from zero with increasing v n 
until a second critical velocity v C 2- There the helix be- 
comes unstable, and the resulting state depends on how 
the flow is applied. A continuously sustained dissipative 
1 texture is found when a constant current is driven in a 
channel. This case has been studied by several groups, 
most extensively by Bozler and collaborators ]IT|-{i"3|] . In 
our rotating cylinder, a stationary state is restored after 
one or more vortex lines are formed in a phase slip. 

Experiment. — Our sample container is a cylinder which 
is rotated around its axis with angular velocity SI. The 
radius of the cylinder R ~ 2 mm is large compared to 
£d, which means that the normal velocity v n — SI x r 
is rectilinear to a good approximation near the cylin- 
drical walls. Therefore, one should see transitions at 
^ci = v c i/R and Sl C 2 = v C 2/R corresponding to the crit- 
ical values of one-dimensional flow. At Sl C 2 dissipation 
sets in only temporarily when a vortex line is created. 
It is driven by the Magnus force into a vortex bundle in 
the center of the cylinder. As a consequence, the effec- 
tive driving velocity v = [SI — Sl v (N)]R is reduced to a 
subcritical value. Here N is the number of vortices in 
the bundle, Sl v = kN/2itR 2 , and k the circulation of one 
vortex line. Compared to channel flow, our experiment 
has high resolution in the measurement of v C 2, since the 
vortices can be counted with a precision of ±5 lines from 
the continuous-wave NMR spectrum. Our sensitivity to 
a helical texture is poor and v c \ has not been observed. 

Three different sample cylinders have been used, which 
were fabricated from epoxy or fused quartz with radii 
R = 2 - 2.5 mm, heights L = 6 - 7 mm, and different 
surface roughness || . No systematic dependence of v c on 
the container was found. A small orifice in the center of 
the bottom plate of the cylinder provides thermal contact 
via a liquid 3 He column to the refrigerator. The NMR 
field H, which is large compared to the dipole field = 
\J 5d/gh ~ 3 mT, is either axial (|| SI) or transverse (_L 
SI). The NMR absorption spectrum has two peaks. The 
frequency shift of the main peak from the Larmor value 
is used for thermometry. The satellite peak arises from 
vortex lines |l4[ and its intensity is proportional to the 
number of lines N. 

Results. — Fig. [I] shows two measurements of the am- 
plitude of the satellite peak as a function of SI. The ac- 
celeration is started from rest (SI = N = 0) at a slow rate 
(dSl/dt = 10~ 3 - 10~ 5 rad/s 2 ), and the satellite intensity 
remains zero until a critical velocity v c = Sl c R. In the 
top frame the amplitude starts to increase linearly when 
fL is exceeded. This means that vortices are nucleated 




FIG. 1. NMR absorption at the vortex satellite peak dur- 
ing a slow acceleration - deceleration cycle: (a) regular vor- 
tex formation at approximately constant v c and (b) burst-like 
vortex formation followed by the regular process. The dashed 
line shows the extrapolation to N — > in the final state where 
the regular process has set in with a reduced ir- 



regularly at a v c that is approximately independent of 
the number of vortices N in the center of the cylinder. 

The acceleration in Fig. la is stopped at a velocity 
!l mM . In order to determine N and v c in this state, the 
rotation is decelerated until at f2 m i„(iV) vortices start to 
annihilate Jr|. The regular parallelogram- like shape of 
the acceleration - deceleration loop in Fig. [j]a shows that 
v c is approximately constant. 

On repeating the measurement, we generally find a 
considerable spread in v c . The distribution is evident in 
Fig. |[ which shows v c measured for different histories 
of sample preparation. Because this variation is much 
larger than seen in one measuring run (Fig. la), it has to 
arise from different metastable states of the system. The 
only source of metastability in our system are different 
superfluid states, i.e. different textures of 1, m, n, and d. 
Additional evidence for the textural origin of the spread 
is discussed below. 

In Fig. 0b the response to acceleration is a sudden jump 
in the signal, which corresponds to a burst of AN w 90 
vortex lines. The critical velocity v c after the jump is 
changed and is generally smaller than before. The sub- 
sequent deceleration to the annihilation threshold in Fig. 
lb shows that v c is reduced from 1.2 mm/s to 0.26 mm/s. 
Obviously, this behavior is caused by a transition from 
one texture to another. The vortex bursts (o in Fig. 2) 
take place only at high temperatures T > 0.7 T c . This is 
consistent with the fact that the energy barriers between 
different textures are smaller at higher temperatures. 

Curve v C 2 in Fig. is a theoretical result for the vor- 
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FIG. 2. Critical velocities as a function of the temper- 
ature T in axial field (H = 9.9 - 15.8 mT, pressures 29.3 
- 34.2 bar). The experimental results, which represent a 
variety of sample histories, are classified as regular (•) and 
burst-like (o) vortex formation. The solid lines denote insta- 
bilities towards vortex formation in different one-dimensional 
textures: v C 2 of uniform or helical texture, u c ls of a locked 
soliton texture, and « c us of an unlocked soliton. The dashed 
lines show the phase boundary between uniform and helical 
textures v c i, and the instability to singular vortices v c sv (es- 
timated from Ref. [[li)). Inset: Distribution of the number of 
vortices AN in a burst, normalized to the maximum possible 
increase AiV max = (2tvR 2 / '/t)(fi b - 0.16-y/ft^) - Ni, where N 
is the initial vortex number and S7b is defined in Fig. |l|. 




FIG. 3. Critical velocities in selected cases, (i) We cool 
slowly through T c under rotation. In axial field this gives re- 
sults (•, ■) near the instability limit n C 2 of the circular texture 
(inset a). In transverse field the corresponding results (o, □) 
are near the instability limit « c y of the expected texture 
(inset c). There is no crucial difference if the magnetic field 
H is on during cooling (■, □) or it is switched on later (•, o). 
The scatter in the data (which exceeds that of v c in Fig. la) 
indicates that the control over the texture is limited, possibly 
due to the end plates of the cylinder, (ii) In transverse field 
v c seems to be less sensitive to the texture, as shown by data 
(A) with varied prehistories. (iii) Often a transverse soliton 
plane (_L Si) is formed when cooling rapidly through T c at 
Q, = 0. In these cases a low v c is found (x, in axial field). 



tex instability of an initially uniform texture. It roughly 
agrees with the largest measured u c 's. Some measured 
values are larger than the theoretical upper limit, which 
may arise from inaccurate parameter values in the cal- 
culation (see below). In order to justify theoretically the 
texture dependence of v c , we have calculated two sim- 
ple cases of initially inhomogeneous texture. A locked 
soliton (LS) is a planar object where both d and 1 turn 
from parallel to v on one side to antiparallel to v on 
the other side of the wall, while they are locked to each 
other (d = 1). Such a texture becomes unstable against 
vortex formation at the velocity i> c ls- Another type of 
domain wall is an unlocked soliton (US), where d = 1 
on one side and d = — 1 on the other. When its plane 
is perpendicular to v, vortices are created at the veloc- 
ity ?; c tjs [0 ■ These simple cases span the spread of the 
measured v c in Fig. 2. This makes it plausible that the 
three-dimensional texture and its defects in the cylinder 
are responsible for all of the variation in v c . 

The data in Fig. || has been collected under a variety of 
prehistories of sample preparation. We now demonstrate 
that more reproducible results are obtained if the initial 
state is prepared using a procedure which favors an equi- 
librium texture. Experimentally textural metastability 
is best avoided by cooling slowly (dT/dt ~ — 1 /iK/min) 
through T c at nonzero £1. For consistent v c values, 



must not be reduced to zero at any point during the v c 
measurements. The critical velocities measured under 
these conditions are shown in Fig. ^ in axial (•, ■) and 
in transverse field (o, □). 

Insets (a) and (b) in Fig. || display two candidates for 
the 1 texture in a rotating cylinder in axial field. The 
"circular" texture (a) is symmetric in rotations around 
the cylinder axis but the "double-half" (b) has only re- 
flection symmetry. In both textures the boundary condi- 
tion fixes 1 perpendicular to the wall but the counterflow 
bends it azimuthal already at distances ~ £d- Concerning 
the flow properties, the circular texture is equivalent to 
the uniform texture, and thus vortices are here expected 
to be created at v C 2- In contrast, the double-half texture 
contains two locked solitons, which means that vortices 
are created already at w c ls- The fact that the measured 
v c data in Fig. || lies consistently above w c ls suggests 
that the circular texture is created in cooling through 
T c in rotation. Thus the measurement of v c gives infor- 
mation about textures which is difficult to get by other 
techniques [ [L8|| . 

In transverse field the texture depicted in inset (c) is 
expected. The flow is here parallel to the field in two sec- 
tors of the cylinder. There it becomes unstable toward 
vortex formation at u c || = u<iP|| / \fp±Xfi±- — P|j)> as calcu- 
lated by Fetter (To). The measured data (o, □, A) agrees 
with this prediction. 
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Fig. H also displays v c measurements (x) in the pres- 
ence of an unlocked soliton, whose plane is perpendicular 
to SI. It is identified from a characteristic satellite peak 
in the NMR spectrum JlSfl , Also this object seems to 
reduce v c to a low but well defined value. 

The measurements of v c can also be carried out as 
a function of the magnetic field. The procedure is to 
accelerate in the desired field, after which the field is 
changed to the NMR value at constant f2, and then v c is 
determined by decelerating to the annihilation threshold 
(f2 m ; n in Fig. [I]). The method works if v c (H) does not 
increase with decreasing H. Indeed we find (axial field, 
33.9 bar, T/T c = 0.75 - 0.85) that v c starts to decrease 
below about 3 mT, but remains finite down to zero field: 
v c (H = 0) = 0.2 - 0.6 mm/s. This is somewhat lower 
than deduced by Saundry et al. from their measured crit- 
ical flow rate at temperatures T/T c > 0.9 

Calculations. — The critical velocities were calculated 
assuming that the order parameter depends only on one 
coordinate x. The energy functional (]l|) was minimized 
numerically for different values of the drive velocity v n . 
The values of the parameters are from Ref. p9[ , except 
that we use determined from the NMR shift in the 
B phase, corrected by trivial strong-coupling effects. No 
adjustable parameters are contained in the calculations. 

For helical textures the minimization gives the energy 
F(v n ,q), where q is the fundamental wave vector of the 
helix. Although the periodic boundary conditions in the 
rotating container tend to fix q, our simulation of the dy- 
namics gives transitions from one value of q to another. 
As a consequence, q approximates qo(v n ), which corre- 
sponds to the minimum of F(v n ,q). (For stability both 
the eigenvalues of the second derivative matrix of F(v n , q) 
have to be positive.) At v C 2 the opening angle (3(x) lo- 
cally starts to grow larger than (3 < 60° found for stable 
helixes and eventually sweeps through 180°. Previously 
v c i and an estimate of v C 2 have been calculated near T c 
by Lin-Liu et al. f|. We find that in high field (Fig. 2) 
the helical texture is not stable above 0.8T C , but at lower 
temperatures v C 2 grows considerably above v c \. 

In low fields (H < H^) u c i drops and vanishes at H = 
when T < 0.85 T c , while v C 2 is nearly independent of 
H. The observed reduction of the measured v c probably 
arises because the texture at small H becomes more sus- 
ceptible to different perturbations such as heat flows or 
end plates of the cylinder. 

In the presence of solitons, the precession of the whole 
soliton texture around x leads to phase slippage. In zero 
field, both w C LS and v c us vanish because there is nothing 
to resist the precession. Such precessing states of both 
the LS and the US have been studied in Refs. p0| , ^l| . 
In contrast, a field H perpendicular to x gives rise to 
a finite v c , as calculated approximately for the US by 
Vollhardt and Maki |2^,|5|]. We find a much lower u c us 
than reported previously. 

Conclusions. — Our measurements of vortex formation 



in 3 He-A are the first to allow detailed comparison with 
theoretical calculations. The quantitative agreement is 
much better than found for v c in other superfluids ( 4 He- 
II and 3 He-B). We find that the critical velocity depends 
on the bulk l(r) texture. The maximal critical velocity 
we associate with the equilibrium texture in axial field, 
while the minimal velocity is a characteristic of textures 
incorporating unlocked solitons. 
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